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Abstract 



The 7] production in the (n, n') bottomonium transitions T(rt) — > 
T(n')r], is studied in the method used before for dipion heavy quarko- 
nia transitions. The widths r„(n, n') are calculated without fitting 
parameters for n = 2,3,4, 5, ml = 1. Resulting 1^(4,1) is found to be 
large in agreement with recent data. 

1 Introduction 

The 7] and 7r° production in heavy quarkonia transitions is attracting atten- 
tion of experimentalists for a long time pp. The first result refers to the 
ip(2S) — * J/ip(lS)rj process (to be denoted as ip(2, l)rj in what follows, sim- 



ilarly for T) with ^ = (3.09 ± 0.08)% [T], T tot = 337 ± 13 keV. 

For the T(2, l)rj and T(3, 1)77 transitions only upper limits B < 2 ■ 10 -3 
and B < 2.2 ■ 10~ 3 were obtained in [2] and [3] correspondingly and prelimi- 
nary results appeared recently in pQ, B(T(2, 1)77) = (2.5±0.7±0.5)10~ 4 and 
B(T(2, 1)tt ) < 2.1 • 10" 4 (90% c.l.). On theoretical side in [5] small ratios of 



widths 



r(T(2,l)fy) 

r(V(2,i)7?) 



2.5 • 10" 3 and 



r(T(3,l)q) 
r(V(2,l)r/) 



1.3-10 



-3 



(1) 
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have been predicted, with the model property that the bottomonium yields 
of r) would be smaller than those of charmonium;specifically in the method 

of [6], the width ratio is proportional to O f(^) J ~ 0.1, for a discussion 

see also [6j [7] . 

However recently [8] new BaBar data have been published on T(4, l)r] 
with the branching ratio 

£(T(4, l)rj) = (1.96 ± 0.06 ± 0.09)10" 4 (2) 

and 

IYT(4 l)n) 

-= 2.41 ±0.40 ±0.12. (3) 



T(T(4, 1)7T + 7T- 

This latter result is very large, indeed the corresponding ratio for ip(2, 1)77 
transition is ~ 0.2 and theoretical estimates (CQ) from [5J for a similar ratio 
of T(3, 1)77/ 7T7T yield 0.015. All this suggests that another mechanism can 
be at work in single r\ production and below we exploit the approach based 
on the Field Correlator Method (FCM) [9] recently applied to T(n,n')7nr 
transitions with n < 3 in [9j [10], n < 4 in [11] and n = 5 in [121 [13] . 

The method essentially expoits the mechanism of Internal Loop Radiation 
(ILR) with light quark loop inside heavy quarkonium and has two fundamen- 
tal parameters - mass vertices in chiral light quark pair qq creation M br w /„. 
and pair creation vertex without pseudoscalars M w « 2u, where uj{uj s ) is the 
average energy of the light (strange) quark in the B(B S ) meson. Those are 
calculated with relativistic Hamiltonian [13] and considered as fixed for all 
types of transitions uj = 0.587 GeV, uj s = 0.639 GeV. 

Any process of heavy quarkonium transition with emission of any number 
of Nambu-Goldstone (NG) mesons is considered in ILR as proceeding via 
intermediate states of BB,BB* ± c.c.,B s B s etc. (or equivalently DD etc.) 
with NG mesons emitted at vertices. 

For one r) or 7r° emission one has diagrams shown in Fig.l, where dashed 
line is for the NG meson. As shown in [HI QUI fTT] . based on the chiral 
Lagrangian derived in [15], the meson emission vertex has the structure 



u = exp y^b—^-J , VaK = V% 



Ccdl = -i I <Fx<iP{x)M br U{x)4>{x) (4) 

(5) 
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(a) 



(b) 



Fig.l Single eta production (dashed line) from Y{n)BB* vertex (a), and BB*T(n') 
vertex (b). 

The lines (1,2,3) in the U matrix ([2]) refer to u,d,s quarks and hence to the 
channels B + B~ , B°B°, B®B® (and to the corresponding channels with B* instead 
of B). Therefore the emission of a single r] in heavy quarkonia transitions requires 
the flavour SU(3) violation and resides in our approach in the difference of channel 
contribution BB* and B S B*, while the 7r° emission is due the difference of B°B°* 
and B + B~* channels (with B — > D for charmonia). 

The paper is devoted to the explicit calculation of single rj emission widths in 
bottomonium T(n, l)rj transitions with n = 2,3,4,5. Since theory has no fitting 
parameters (the only ones, M u and M& r are fixed by dipion transitions) our pre- 
dictions depend only on the overlap matrix elements, containing wave functions 
of Y(nS), B,B S , B* , B*. The latter have been computed previously in relativistic 
Hamiltonian technic in p3] and used extensively in dipion transitions in jll|, ll2 ^ [T3] . 

The paper is organized as follows. In section 2 general expressions for process 
amplitudes are given; in section 3 results of calculations are presented and discussed 
and a short summary and prospectives are given. 



The process of single NG boson emission in bottomonium transition is described 
by two diagrams depicted in Fig.l, (a) and (b) which can be written according to 
the general formalism of FCM [9j [TT| [T2] as (we consider rj emission) 



2 General formalism 



M = + M®\Mf = M% Bt - M 



BB*i 



i = 1,2 



(6) 
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(2) 

where has the same form, but without NG boson energy in the denominator 
of ©. Here 7 = , M " M '"- 2 and the overlap integral of T{nS) and £5* wave 

functions is (for details see Appendix) 

jW(p,k) = 4 3 (°)/„,^(p)e" fc "% (8) 



Jn'(P) = ^2 (1) VfiB.(p)e ^ (9) 

z„(BB*) and Z2{BB*) are Dirac traces of decay matrix elements T(nS) — > BB*r/ 
and — > T(n'), respectively they are defined in [12 [H] and below in Appendix. 
The special point in our case is that r\ meson is emitted in P wave, hence one must 
extract the corresponding term in the Dirac trace, for details see Appendix. 

z\BB ).z 2 (BB ) = i^—j — g^s— A*. (10) 

Here u n = A ~ sJt> an d ^(^s) is the average energy of the b quark 

in B{B S ) meson; from Table IV in [9] one finds that Q = 4.827 GeV, Q s = 4.830 
GeV. In what follows we shall neglect the difference between 0, Q s and the mass of 
b quark = 4.8 GeV. Note that these large masses cancel in all matrix elements 
and final expressions will depend only on energies u and lo s and differences of 
threshold positions: AM* = M(B) + M(B*) - M(T(nS)) and AM S - the same 
for B s , B* masses. Note, that the contribution of the B*B* , B*B* channels vanish 
hence we shall consider only BB* and B S B* channels. 

Indices i'i in e^n in ()10p refer to the T(n'S') and T(nS) polarizations respec- 
tively. Finally, coefficients faifli and A n = 201 + refer to the expansion of 
realistic wave functions of T(nS), T(n'S) and B, B* , B s , B* s computed in [T3] in 
series of oscillator functions and (5\ , P[ and denote the \ 2 fitted oscillator pa- 
rameters for those functions respectively, see [11] for details. 

Finally we define all quantities in the denominator of ([7]); in M BB » the de- 
nominator is 

E-E(p) = M(T(nS))-(u; v +M B +M* B +^-+ iP 2 ~^ 2 ) = -AM*- Wl) -£(p,k). 

(11) 

For Aiii one omits u v and k in (jlip . Finally one can represent the matrix element 
M$ as follows: 

M (i) _ ~ P ..,, k ,jL ( —rW - —tfX^ do) 
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3A„, s 



£(2)__L £ (2) e 4 A 



(13) 



with 



d 3 pp 2 (°)j w (p) W/^e 
( 2 ^ 3 (AM*+n + p2 I (P' k ) 2 



p 2 



r (2) = /" <*W Win(g) ^# ^ 
(2.)3 ( AM * + ^ 



(14) 



(15) 



For £i 1} , Cf ] one replaces AM* with AM* and M B , M B with M Bs , Mb* . Here 

^0 - An + A„; • 

The width of the Y(n,n')r) decay is obtained from \Ai\ 2 averaging over vector 
polarizations as 

(6P £«\ (cf ] 1 



1 OP -— n 



27 



U} 3 LU 3 



+ Un' 



U 3 UJ 3 



2 (16) 

where the phase space factor d<& = -^j2irS(M(T(n))-M(T(n'))-ui v - 2 M(T(n')) )• 



Introducing the average w = + w), one can rewrite (fl~6|) as 



M &r \ /M 



2 - — 2 j.3 
2lu J LO^ 



u„> 







K 








a) 









£(1) 



+ 



(17) 



with <, = ^ 7 • lO" 8 . 



One can see from the general structure of r„, that the main effect comes from 

:; -s |0.882 - 1.139| « 0.257, and from the difference of 



the difference 



|£« <o.05. 



3 Results and discussion 

We consider here the single rj emission in bottomonium transitions T(n, l)rj with 
n = 2, 3, 4, 5. The corresponding values of AM*, AM*, uj„, k are given in the Table 
1. 
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Table 1. 

Mass parameters of T(n,n')rj transitions (all in GeV, k in GeV/c) and matrix 
elements C^,cf (in GeV). 



(n, n') 


2,1 


3,1 


4,1 


5,1 


AM* 


0.582 


0.25 


0.026 


-0.26 


AM* 


0.757 


0.425 


0.20 


-0.08 




0.562 


0.87 


1.075 


1.325 


k 


0.115 


0.674 


0.923 


1.20 


£(1) 


0.263 


-0.188 


0.29 


1.48 -10~ 3 


& 


0.240 


-0.174 


0.121 


1.41 -10" 3 


£(2) 


0.390 


-0.340 


0.255 


-0.347 
-i0.0785 


/•(a) 


0.341 


-0.298 


0.226 


+0.0584 
+i0.102 



The resulting values of r^(n, n') have been computed as in (fTTl) with to = 0.587 

GeV and uj s = 0.639 GeV, calculated earlier in [14], see Table 4 of [9], and u n = 
with /?2 = 0.48 GeV, and A n both fitted to the realistic wave functions in [11], with 
A n = 2.56; 1.54, 1.21, 1.05 and 1.35 (all in GeV 2 ) for n = 1,2,3,4,5 respectively. 

Results of calculations are given in Table 2. 

Table 2. 

Values of T v (n,n') (in keV) calculated using Eq. (fTTl) vs experimental data 
T^{n,n') (in keV). 



(n, n') 


2,1 


3,1 


4,1 


5,1 


r„(n,n') 


5.0-10" 2 


2.9 


1.81 


7.04 






(0.8±0.3) • 10~ 2 
[4] 




4.02 ± 0.6 

[8] 





Looking at the Table 2, one can see, that there is an order of magnitude 
agreement with experiment. Indeed, the factor ( ; y i21 ) flif) can ^e estimated 
from T(n, n^nir transitions studied in [5]-|13] to be roughly in the range [|, 2]. At 
the same time r^ xp (2, 1) differs from our calculated value several times, and more 
accurate measurements as well as theoretical calculations are highly welcome here. 
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Another point, not shown in Table 2, is the old upper limit [3] on -6^(3, 1), namely 
5(1(3,1)77) < 2.2- 10" 3 [3], which yields 1^(3,1) < (4.5 ±0.4) • 10" 2 keV and is two 
orders of magnitude below our calculated value. Hopefuly new measurements can 
resolve this disagreement. On theoretical side our formulas (|14p - (|16p automatically 
produce the width T r? (n, n') of the order of 0(1 keV), for all (n, 1) transitions 
except for (2,1), where a small phase space factor k 3 gives two orders of magnitude 
suppression of T r? (2, 1). For the T ri (5, 1) one obtains a 7 keV value, which is however 
small as compared with the r 7rvr (5, 1), the latter being 0(1 MeV). For T r? (4, 1) and 

r 7r7r (4, 1) from p] the calculated ratio is R v/n7T = = 3 (Jfc) ~ 3 

which roughly agrees with experimental value R C ^ T7T = 2.41 ± 0.40 db 0.12. 

To check stability of our results, we have used for the wave function of B s 
the realistic wave function different from that of B. As a result one obtains for 
-. V \ {n i ' l) the values (2.74 • 10" 2 ; 1.13; 0.44; 7.3) keV for n = 2,3,4,5 respec- 

( Hbr\ ( M " ^ 2 
\ fir J \ 2^ ) 

tively, which should be compared with numbers in the upper line of Table 2. The 
same type of sensitivity occurs for modifications of other wave functions, imply- 
ing that our results strictly speaking yield the correct order of magnitude but not 
exact values of Y v . 

Summarizing, we have calculated the single 77 production width r r) (n,n / ) for 
T(n, 1)77 transitions with n = 2,3,4,5. We have found that T r] (n, 1) are of the 
order of and larger than T n7T (n, 1) for n = 3,4. This fact is in agreement with 
the latest measurements in [8] of T^ xp (4, 1) and disagrees with earlier experimental 
limit on r^ xp (3, 1). Our calculations do not contain fitting parameters; the only 
two parameters M& r , M u are fixed by previous comparison with dipion data. One 
should stress that 77 production in bottomonium is not suppressed in our approach 
as compared to 77 production in charmonium transitions. This is in contrast with 
the results of method of [6]. 

The financial support of grants RFFI 06-02-17012, 06-02-17120 and NSh- 
4961.2008.2 is gratefully acknowledged. 
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Appendix 1 

Matrix element of single r/ emission 

According to the general theory in [91 [11], the matrix elements M.^\m.^ 
for T(n, n')rj corresponding to diagrams of Fig.l, (a) and (b) respectively, can be 
written as 

(1) MbrMu, f d 3 p v Jnn^Un'n 2 n 3 (P) ^ 

^ 1 > f v N c ^J-J (2vr)3 n ^ 3 E-E n2n3 (p) ' 1 ' ) 

Here n-2, 113 are channels of intermediate state, with e.g. ri2 = B, B* , B s , B*, we 
omit indices 712 , re 3 and write 

J n (p, k) = J AS_§ n ( cp - I + q)V5„ 2 (q)^n 3 (q - k) (A.2) 
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where *ff n ,ifi ni are momentum space wave functions of T(nS) and B{B*) mesons 
respectively. 

The vertex factor y ^3 is calculated in the same way as in [9] , namely from the 
Dirac trace of the projection operators for the decay process, in our case this is 
T(nS) — > BB*rj. Identifying the creation operators as "^bli^bi ^ftTs^Vu tyblk^m n = 
it, d, s and extracting vertex of 77 creation from the Lagrangian AC = — J ^rJJ r Mi, r '^ n d 4 x 

(I \ 

which gives i—trls^Mit^ with A = -4= 1 , one has for the decay pro- 



fn y 2^77 

cess (cf. Appendix 1 of [9]) 

G(T -> 55*7/) = tr[jiS b (u, w)j 5 Sn(w, x)j 5 S n (x, w)^ k S- b (w , v)\ (A. 3) 



-2 



As shown in [9], appendix 1 and 2, the quark Green's functions can be split into 

two factors S(x, y) = A ± G(x, y), with the projection operators A^ = mk — 

k = b,n and the scalar part G(x, y), where spins are present only in spin-depdndent 
interaction and treated as corrections. Here lo^ is the average energy of quark in 
given meson. Hence one is brought to the spin factor Z. 



Z = *r( 7i A+75A-7 5 A+7 fe A 6 -) (A.4) 

which is equal to 

z = ^±j£ ((^ 2 - pV)** - pH + pIp!)- (A.5) 

Here r2,w are average energies of 6 and n quark in B or B* . One can identify 
the momenta of B and B* as Pi = p and P2 = — p — k, then q in (|A.2[) can be 
expressed as 

Pg = -q+ — -xP, P q = q T7TP~ k lT^ ' ( A6 ) 

W + li H UJ + \l \L + UJ 



and Z is (we put mj, = £1) 



Z = ^ (-kq5j fe + fc^ + fc fc g;). (A.7) 



It is important, that we are looking for the P-wave of emitted ij, and hence for P 
wave of relative BB* motion, hence the integral (|A.2p should yield the term pk. 
This indeed happens, when one approximates ^ n ,ip n as series of oscillator wave 
functions and (|A.2|) has the form 

PL k 2 

J n ( P ,k) = ^ 23 e A " ^ (°)/ n (p). (A.8) 
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In the process of dq integration in ()A.2j) one changes the integration variable 
qi — > q[ — u n pi + O(ki) with u n = /3|/A„ are oscillator parameters, found by x 2 
procedure. 

Thus result of d 3 q integration yields 

2/123 = ^(- k P^fc - hPk + hPi)- ( A -9) 
In an analogous way one obtains for J n '(p) in (jA.ip the form 

MP) = &e~ p2/A "' W Mp) (A.10) 
and y^}\ 3 is obtained from the Dirac trace for the process BB* — > T(n ; S), 

Z 2 {BB*) = J- em (-2 qi + -^Pi) (A.11) 

and the result of integration over d?q yields in (jA.lOj) 

& = -ei'ki-- (A-12) 

UJ 

Here i 1 is the polarization of T(n'S) (represented by ipbli'^b) an d k as in (|A.9P 
is the polarization of B* . Averaging over angles of p one obtains 

(&y%) P = Y^ e ^ (A.i3) 

and finally one writes as in (|12p 



For A4^(n, n ) one can use time inversion and interchange indices and 
change sign of k, obtaining in this way Eqs. (I14p and (I15p of the main text. For 
the intermediate state of B*B* the summation over polarizations of B* yields a 
net zero result, therefore we are left with only (BB* + B*B) intermediate state. 
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